MZ-TH/10-45, EFI Preprint 10-30 



Long-Distance Dominance of the CP Asymmetry in B — y X Sj d 7 Decays 
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We show that in the Standard Model the parametrically leading (by a factor l/a s ) contribution to 
the inclusive CP asymmetry in B — > X s ,d 7 decays arises from a long-distance effect in the interfer- 
ence of the electro-magnetic dipole amplitude with the amplitude for an up-quark penguin transition 
accompanied by soft gluon emission. Using model estimates for the associated hadronic parameter 
Ai 7 , we predict a value in the range —0.6% < •Ajcf-, < 2.8%. In view of current experimental data, 
a future precision measurement of the flavor- averaged CP asymmetry would signal the presence of 
new physics only if a value below —2% was found. A cleaner probe of new physics is offered by the 
difference of the CP asymmetries in charged versus neutral B — y X s y decays. 



The radiative decays B — y X s j and B —y Xjy are im- 
portant for testing the Standard Model (SM) and probing 
its possible extensions. Both the CP-averaged branching 
ratios and the CP asymmetries are useful in this con- 
text. Thanks to a vigorous effort, the experimental error 
on the B — > X s ^f branching ratio has been reduced below 
10%, which is close to the (irreducible) theoretical uncer- 
tainty in the calculation of this observable. For the CP 
asymmetry, on the other hand, the experimental error is 
thought to be much larger than the theoretical one. The 
current world average is [l| 
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T(B -» X sl ) - T(B -y X- sl ) 
Y{B -y X sl ) + T(B -y X s j) 



-(1.2 ±2.8)%. 



(1) 

It is widely believed that the theoretical prediction for 
the CP asymmetry in the SM is short-distance domi- 
nated, leading to a tiny value of about 0.5% as a re- 
sult of a combination of perturbative, CKM, and GIM 
suppression 0-0] • A dedicated analysis finds A^f^ = 
(0.44 ± °;^±0.03 1 ^)% @, where the errors refer to un- 
certainties associated with the quark-mass ratio m c /mb, 
CKM parameters, and higher-order perturbative correc- 
tions. This suggests that finding an asymmetry outside 
the range < Ax,-y < 1% would be a clean signal of 
new physics. Indeed, various extensions of the SM have 
been analyzed with regard to the constraints ari sing from 
the measured value of the CP asymmetry 0, 0- 15 1 , and 
reducing the experimental error to the 1% level is consid- 
ered an important goal of future super-flavor factories. 

In recent work [l6j , we have presented a new factoriza- 
tion formula for inclusive radiative B decays in the rel- 
evant region of large photon energy. In addition to the 
familiar "direct photon contributions" , in which the pho- 
ton couples to a local operator mediating the weak decay 
in the effective low-energy theory, novel "resolved photon 
contributions" appear. They account for the hadronic 
substructure of the photon, which is probed when it cou- 
ples to light collinear partons. Importantly, even after 
integrating over the photon energy spectrum, the re- 



solved photon contributions cannot be described using 
a local operator-product expansion. They give rise to 
first-order Aqcd corrections to the inclusive decay 
rate and CP asymmetry, which are proportional to B- 
meson matrix elements of non-local operators in heavy- 
quark effective theory (HQET). The analysis of 16| has 
shown that these non-perturbative effects lead to an irre- 
ducible uncertainty in the theoretical prediction for the 
CP-averaged B — y X s j branching ratio of about ±5%. 

An interesting feature of the resolved photon contri- 
butions is that they give rise to novel, calculable strong- 
interaction phases related to hard-collinear jet functions, 
which are convoluted with real, non-perturbative soft 
functions. It is interesting to explore the potential impact 
of these contributions on the inclusive CP asymmetries. 
While these effects are still of order Aqcd /nrib, we show 
that they give the parametrically leading contribution to 
the B — y X s ^d 7 asymmetries in the SM, and that they 
can strongly influence the way in which new-physics ef- 
fects might show up. 

Experiments measure the CP asymmetry Ax s ^{Eq) 
defined with a lower cut on the photon energy, _E 7 > Eq, 
with Eo ranging between 1.9 and 2.2 GeV. Detailed the- 
oretical studies have shown that the dependence of the 
asymmetry on the value of Eq is very mild [3|. We will 
assume for simplicity that the cutoff can be chosen suf- 
ficiently low, such that A = — 2Eq = few x Aqcd is 
in the perturbative domain. We will refer to the asym- 
metry defined with such a cut as "partially inclusive" . 
In this case the direct photon contributions to the CP 
asymmetry can be calculated in terms of local operator 
matrix elements usin g a combined expansion in powers of 
A / mi, and Aqcd / A |l7[ . The resolved photon contribu- 
tions can still not be expressed in terms of local matrix 
elements. However, the relevant soft functions can be 
simplified in this limit, as described in [l6j . 

The direct photon contributions to the partially inclu- 
sive CP asymmetry reduce to perturbative expressions 
depending on a cut parameter S = A /m&. At first non- 
trivial order in a s , one obtains 0, [HI 
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where z = (m c /m b ) 2 , and e s = (V ub V* s ) / (V tb V t * s ) = 
\ 2 (if)~p)/[l — \ 2 (l-p + if])] + 0(\ 6 ) in terms of Wolfen- 
stein parameters. The Wilson coefficients of the electro- 
magnetic and chromo-magnetic dipole operators in the 
effective weak Hamiltonian are denoted by CV 7 and Cs g , 
while C\ is the coefficient of the dominant current-current 
operator. Contributions from other operators are negli- 
gibly small in the SM. We suppress the scale dependence 
of a s and Cj. Explicit expressions for the functions v(z), 
b(z,S), and b(z,6) can be found in @, In the SM 
the Wilson coefficients are real, and only terms in @ 
proportional to the imaginary part of e s contribute. The 
numerically most important contributions arise from the 
interference of charm- and up-quark penguin graphs with 
virtual or real gluon emission (first two diagrams in Fig- 
ure [T|) with the leading electro-magnetic dipole ampli- 
tude. The above expression is however more general, as 
it holds for all new-physics models in which the dominant 
non-standard effects are described by additional (possi- 
bly complex) contributions to the dipole coefficients C77 
and Cs g - It would be straightforward to extend the anal- 
ysis to include opposite-chirality dipole operators [3| . 

The above expression simplifies considerably if we 
adopt the power counting m 2 — O (m b Aqcd) for the 
charm-quark mass and expand in powers of z,<5 = 
£^(Aqcd/w6). The terms proportional to b(z,S) and 
b(z,S) scale as (Aqcd/w;,) 2 and can be neglected to a 
good approximation, while the contribution proportional 
to v(z) contributes at first order in Aqcd/™& and can be 
simplified by expanding v{z) to zeroth order in z. This 
yields the approximation 



,idir 



— Im— 

81 C77 



9 Im cv; 



(3) 



40A C 
9m b 



Im 



(1 



cv 7 




which is independent of the cutoff Eq. Here we have 
introduced the scale A c ~ Aqcd defined as (we use m b — 
4.65 GeV and m c = 1.13 GeV as in (HI) 
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0.38 GeV. 
(4) 

In the SM only the last term in ((3]) contributes, which 
exhibits the triple suppression by a s , Im(e s ) ~ A 2 , and 
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FIG. 1: Penguin diagrams with virtual and real gluon emis- 
sions. The gluons in the first two graphs are highly energetic. 

(m c /m b ) 2 ~ Aqcd/^6- 

Using the factorization analysis of B — > X s j decay per- 
formed in 16], it is possible to derive for the first time 
the resolved photon contributions to the partially inclu- 
sive CP asymmetry. They arise from the interference of 
the electro-magnetic dipole amplitude with resolved pho- 
ton contributions involving up- and charm-quark penguin 
transitions or chromo-magnetic dipole transitions. The 
result can be expressed in terms of three hadronic pa- 
rameters Aij related to convolution integrals over two 



soft functions denoted by hn{ui) and h 7s (wi, 102). 
lowest order in a s and l/m b we find 



At 



jres _ 

Ax* 7 — — 
m b 



Im 



(1 + ^)7^ 



Ci 

C77 





' Ci " 


A c 17 - Im 


e s — 

L-7 7 _ 



A 



17 



l m 4 7rQ , s Af 8 



G 



(5) 



where (omitting the scale dependence of the soft func- 
tions and parameters) 
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The functions hij are defined in terms of S-meson matrix 
elements of non-local operators in HQET jl6j| . In light- 
cone gauge n ■ A — 0, we have 

f dt (Bim^npgG^MhmB) 
h 17 (u>) = J-e — , 

(8) 
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where — (1,0,0,1) and — (1,0,0,-1) are two 
light-like vectors, and h are effective heavy-quark fields 
in HQET [20| , Similarly, /i 7 g (wi, 0J2) is given in terms of 
a matrix element of a four-quark operator (l9j . The func- 
tion hn is an even function of u, whose normalization is 
equal to the HQET parameter 2A 2 « (0.5 GeV) 2 . No rig- 
orous constraints are known for the function h 7 \ . Note 
that the parameter Afg in ([5]) depends on the flavor of the 
spectator quark inside the B meson. In the limit of SU(3) 
flavor symmetry, it can be shown that Afg = e spC c A7g 
fl6j | . where e spe c denotes the electric charge of the spec- 
tator quark in units of e (e sp0 c = 2/3 for B~ and —1/3 for 
B°). Evaluating the hadronic matrix element of the cor- 
responding non-local four-quark operator in the vacuum 
insertion approximation, we find that 



A 7 g ~ e spoc A78 



^spcc 
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where Jb ~ 193 MeV is the B-meson decay constant and 
0+ its leading light-cone distribution amplitude (2lj . 

At present, there does not exist any systematic the- 
oretical approach to determine the hadronic parameters 
Ay from first principles. Numerical estimates must then 
be obtained by modeling the correspon ding soft func- 
tions. Employing the models studied in [16], we obtain 
the ranges 



330 MeV <A? 7 < +525 MeV, 



-9 MeV <A C 17 < +11 MeV, 
17 MeV <A 78 < 190 MeV. 



(10) 



All three estimates are very uncertain, but we observe 
that A" 7 and A 78 are expected to be of order Aqcd • The 
slight preference for positive values of A\ 7 is due to the 
normalization constraint on the function hyj mentioned 
above. Note that in the formal limit m c — » m u = 
the values of A" 7 and A\ 7 coincide. However, we pre- 
dict a strong GIM violation owing to the fact that the 
integral in the second relation in © starts at 4m;? /mj, 
1.1 GeV, at which the soft function is expected to take 
already rather small values, since it is governed by non- 
perturbative dynamics and must vanish for u> — > 00. 

The complete theoretical result for the partially inclu- 
sive CP asymmetry in B — > X s ~/ decay is obtained by 
adding the direct and resolved contributions (J2)) and (|5|) . 
In order to understand better the structure of the result, 
we now present some approximate formulae obtained by 
using expression ([3]) for the direct contribution. Our nu- 
merical results will always be derived using the exact 
expression. For the CP asymmetry in the SM, we obtain 
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+ 0.71 % . 



We fix the photon cut at E = 1.9 GeV and use A = 
0.2254, p = 0.144, f\ = 0.342 for the Wolfenstein param- 
eters. Our choice \x = 2 GeV for the factorization scale 
(for which a s (n) = 0.307, and d(/i) = 1.204, CV 7 G") = 
—0.381, Cs g ([J.) — —0.175 at leading order) is motivated 
by the fact that the strong phases required for a non-zero 
CP asymmetry arise either from GIM violations related 
to charm-quark loops (for which fi ~ 2m c ), or from cut 
hard-collinear propagators (for which fj, ~ \A™&Aqcd). 
In the resolved photon term we keep the contribution of 
A^ 7 in order to make explicit that the CP asymmetry 
vanishes in the formal limit m c — m u due to the GIM 
mechanism. In practice, however, this contribution can 
be safely neglected. The dominant contribution arises 
from the up-quark penguin graph with emission of a soft 
gluon, in which the quark loop is probed at a light-like 
distance away from the heavy quarks (last diagram in 
Figure [lj . We do not show the dependences of our re- 
sults on variations of the input parameters, which are 
much smaller than the uncertainties associated with the 
resolved photon contributions. Our central value 0.71% 
for the direct photon term is larger than that obtained 
in [B[ since we use smaller values for fj, and m c . 

The resolved photon term proportional to A^ 7 in (fTTj) is 
parametrically larger than the direct photon term, which 
contains an additional a s suppression. Numerically, this 
term dominates as long as |A" 7 | is larger than about 
200 MeV. Using the model estimates shown in (|10[) we 
find the range -0.6% < Af^ 7 < 2.8%, which covers most 
of the experimentally allowed range ([lj . Only a value of 
the asymmetry below —2% could be interpreted as a sign 
of new physics, as in this case A" 7 < —700 MeV would be 
much larger in magnitude than our model expectations. 

In extensions of the SM, in which the dipole coefficients 
CV 7 and Cs g receive new CP-violating contributions, ad- 
ditional terms arise. Using the approximation ([3]), we 
find 
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For the first two terms the resolved photon contribu- 
tions give rise to power corrections to the direct photon 
terms, which are numerically significant since a s /-K ~ 
Aqcd I nib- For the third term, which is the only one 
present in the SM, the resolved photon contribution is 
likely to be more important than the direct photon term. 

To illustrate the impact of the resolved photon terms, 
we consider the class of new-physics models in which the 
dominant non-standard effects are encoded in the val- 
ues of the dipole operators, which we parameterize in 
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the form C 7l /Ci = {C 7l / d) SM r 7 and C 8g /Ci = 
(C 8g /Ci) SM r a e ies . Using that arg(-e s ) « -7 in the 
SM, we then obtain for our default parameters 
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For the flavor-averaged CP asymmetry, i (^4 x - 7 +^l x j7) 1 
we must replace e spC c — > ^- Then the resolved pho- 
ton contributions are subdominant except for the second 
term, which is already present in the SM. In principle 
very large asymmetries are possible from the first and 
third terms [3j, which however are already ruled out by 
the data. Once we assume that the effects of new physics 
are at most a few percent in magnitude, it will be diffi- 
cult to disentangle them from the hadronic uncertainty 
due to the A" 7 parameter. 

A non-trivial feature of our analysis is that the resolved 
photon contributions induce a flavor-dependent term in 
the CP asymmetry already at order Aqcd/^6- In the 
SM such effects are suppressed, compared with (|11[) . by 
at least one additional factor of Aqcd/to{, and are thus 
bound to be negligible. We thus propose a future preci- 
sion measurement of the CP asymmetry difference 
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as a sensitive probe for flavor physics beyond the SM. 
Even though it will be difficult to determine the value 
of the hadronic parameter A78 with any reasonable ac- 
curacy, we observe that if either the electro-magnetic or 
the chromo-magnetic dipole coefficients (or both) receive 
a sizable CP-violating new-physics phase, the difference 
AAx s -f can eas iby reach the level of 10% in magnitude. 
It is important in this context that the Wilson coefficient 
of the chromo-magnetic operator can be much enhanced 
with regard to its SM value, so that r%/r 7 ~ a few is 
possible 0. 

All of the above expressions also apply to the CP 
asymmetry in B — > X^j decay, once we replace e s by 
e d = (V ub V: d )/(V tb V t * d ) = (p - ifj)/(l -p+ if,). As 
a result, the CP asymmetry for B — > Xd'y decay in 
the SM differs by that in B — > X s ~f decay by a factor 



Im(€(j)/Im(e s ) ~ —22. With the parameter values shown 
in {10]), we obtain the range -62% < A^ < 14%. It is 
worth emphasizing that in the endpoint region of large 
photon energy there do exist contributions to the photon 
energy spectrum which could affect the CP asymmetries 
in the decays B —> X s j and B —> X^ in a way dif- 
ferent from the above rescaling. They are sensitive to 
four-quark soft functions with flavor content bssb and 
bddb, respectively, whose matrix elements between B d 
meson states can be different even in the flavor SU(3) 
limit. However, as shown in [TJ| these effects integrate 
to zero in the partially inclusive CP asymmetries at or- 
der Aqcd/™6- As a result, the untagged CP asymmetry 
for B — > X s+ dj decays vanishes in the SM (up to tiny 
U-spin breaking corrections @, 0, [22| ) even after the re- 
solved photon terms are taken into account. 
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